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ABSTRACT 



The oscillatory forces generated by the interaction of 
a train of gravity waves with a rigid horizontal circular 
cylinder located near a simulated ocean floor were deter- 
mined experimentally. A circular cylinder of six inch dia- 
meter was subjected to gravity waves of wave lengths varying 
from 2 to 18 feet in a two-dimensional wave channel. Both 
the horizontal and vertical components of the hydrodynamic 
force were measured. 

Experimental results are plotted in dimensionless 
form. The horizontal and vertical dimensionless force co- 
efficients are plotted as a function of the relative water 
depth, the relative wave lengtii, <3.nd tlic relative vvave 
height. It was found that the horizontal force coefficient 
varied quite linearly \\rith the wave height, and the maximum 
value \>ras approximately equal in either direction. The 
vertical force component, however, showed a much more com- 
plex variation. The upward and dovmward magnitudes were 
approximately equal for small amplitude waves, but as the 
wave height increased, the variations became quite non- 
linear and the magnitudes in the two directions differed 
considerably . 
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NOMENCLATURE 



Symbol 



Description 



m 

d 

"x 

X 

+ 

y 

+ 



g 

H 

h 

k 

L 

1 

T 

t 

u 

u 

V 

V 



radius of circular cylinder 

added-mass coefficient, dimensionless 

depth of submergence of cylinder 

maximum horizontal force per unit length 

F / [pga^ (H/2a) ] , horizontal force coefficient, 

^ dimensionless 

maximum positive vertical force per unit length, 
(upward) 

F^/ [,pga^ (H/2a) ] , positive vertical force, coeffi- 
^ cient, dimensionless 

maximum negative vertical force per unit length, 
(downward) 

F / [pga^ (H/2a) ] , negative vertical force coeffi- 
^ cient, dimensionless 

acceleration due to gravity 

height of the incident wave 

water depth 

wave number= 2ir/L 

wave length of the incident wave 

length of the cylinder 

period of the incident wave 

time 

component of fluid velocity in the x-direction 
component of fluid acceleration in the x-direction 
component of fluid velocity in the y-direction 
component of fluid acceleration in the y-direction 
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Symbol 



Description 



X space coordinate 

y space coordinate 

velocity potential 
p density of the fluid 

a frequency of the incident wave= 2 tt/T 

y dynamic fluid viscosity 
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DIMENSIONLESS NOTATION 



d/a 



[Pga^ (H/2a) ] 
Fy/ [Pga^ (H/2a) ] 
Hl2a 
h/a 



Zva/L 



relative depth of submergence 
horizontal wave force coefficient 
vertical wave force coefficient 
relative wave height 
relative water depth 
relative wave length 
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INTRODUCTION 



I . 

In recent years, there has been an accelerated interest 
in the exploration of the ocean bottom with special atten- 
tion devoted to those depths encompassed by the Continental 
Shelf. This interest has been motivated in great part by 
the desire to exploit the mineral assets of the ocean bot- 
tom. Additionally, investigation of ocean bottom pipelines 
as specifically applied to sewage outfall lines, oil pipe- 
lines and large aqueducts are under consideration. 

As a result, there has been considerable interest in 
the hydrodynamic interaction of gravity waves with large 
submarine structures resting on the ocean bottom. The 
design of large submarine offshore structures necessarily 
requires one to have an understanding of the forces induced 
by such environmental phenomena as surface waves. To this 
point, the theoretical work has yet to be a conclusive tool 
in assisting the design engineer to predict all expected 
conditions to which the anchored structure will be subjected 
as it lies on the ocean floor. Experimental observations 
of the forces induced by ocean waves are difficult to ac- 
quire, and it is therefore necessary to simulate as pre- 
cisely as possible the conditions existing at the ocean 
bottom . 

Extensive experimental work has been carried out since 
Stokes, who in one of the earliest papers on the subject 
[1], showed that, by experimentation with a pendulum- type 
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system, the expression for the force on a cylinder oscilla- 
ting in an infinite fluid consists of two terms, one in- 
volving acceleration of the cylinder, and the other 
involving velocity. 

Wave forces exerted on piles have been intensively in- 
vestigated during the past two decades. It can be shown 
that when the relative size of the body, i.e., the cylinder 
diameter, is small compared to wave length, a condition 
which is generally satisfied in the case of piles, the 
familiar Morison equation [2] is a valid approach. This 
equation combines both a drag term and an inertia term, 
and the values of the experimental coefficients of drag 
and inertia are left to the ingenuity of the designer. At 
best, these tv.'o experimental ceefficients arc dependent 
upon the geometry of the objects under consideration and 
the amplitude of the fluid motion, and it would appear that 
20 years of research has shed little light on understanding 
this interdependence. Moreover, the proximity of the free 
surface as well as the impermeable bottom surface affect 
these coefficients. In view of these difficulties^ the 
Morison equation is not relied upon in the present approach 
to the problem, and the results are represented directly in 
dimensionless form rather than through the drag and inertia 
coefficient . 

Keulegan and Carpenter [3] used periodic standing 
waves in their analysis of a cylinder subjected to a field 
of oscillatory motion. They were able to separate the total 



12 



resistance into drag and inertial components. They correla- 
ted values of inertia and drag coefficients under conditions 
of varying current intensity and cylinder sizes with a peri- 
od parameter, U^^^T/D, a dimensionless group combining the 
maximum fluid velocity the wave period T, and the cylin- 

der diameter. This parameter is a form of the relative 
fluid displacement since U^T is proportional to the ampli- 
tude of the fluid motion. 

Sarpkaya and Garrison [4] studied both theoretically 
and experimentally unidirectional flow past a circular cyl- 
inder with constant acceleration. They also showed that 
the total resistance as well as the drag and inertia com- 
ponents of force may be represented as a function of rela- 
tive fluid displacement. Tn the present analysis of the 
experimental data, this same general kind of dependence on 
fluid displacement is supposed, and, accordingly, the force 
coefficients are correlated with the ratio of the wave 
height to cylinder diameter. 

Johnson [5] has investigated the horizontal forces 
resulting from the interaction of gravity waves with a 
circular cylinder placed relatively near the bottom. Since 
only long waves were studied, the resulting forces were 
supposed to be dependent on the water depth and wave height 
only. This study shed little light, however, on the basic 
problem .• 

Herbich and Shank, in a recent study [6], have con- 
ducted wave channel testing of vertical and horizontal loads 
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on submerged three-dimensional structures due to interaction 
with gravity waves. Several models were studied using a 
two-dimensional wave tank, and they concluded that the re- 
sults of these experimental studies were sufficient to be 
applied to preliminary designs. Unfortunately, few actual 
experimental points were presented and therefore the degree 
of correlation is somewhat suspect. 

Garrison, et at. [7], investigated the wave/structure 
interaction problem as applicable to larger objects where 
the Morison equation is no longer valid and a diffraction 
theory must be employed. In this situation, the incident 
wave is scattered upon encountering the large structure, 
and the assumption that the object does not affect the in- 
cident wave is no longer valid. The diffraction theory 
accounts for the relative size of tlie body and, neglecting 
viscous effects, assumes the existence of a velocity poten- 
tial. When large structures, such as the recently deployed 
Chicago Bridge and Iron Persian Gulf oil storage tank, are 
designed, the analysis must include consideration for rela- 
tive size and free surface effect, and this basic consid- 
eration is inherent in the diffraction theory. 

It appears that there will continue to be heavy re- 
liance on experimental model results in the design process 
until more reliable theoretical predictions of wave induced 
forces are assured. To this time, large structures have 
been placed on the ocean floor only after extensive model 
studies have been carried out. 
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The present study was undertaken to obtain an insight 
into the interaction of gravity waves ivith large submarine 
structures located at or near the ocean floor. The study 
includes direct measurement of horizontal and vertical 
forces on the submerged cylinder. The results are pre- 
sented in dimensionless form as a function of the relative 
water depth, the relative wave length, and the relative 
wave height. 

The interest in attempting this investigation was 
generated while participating in numerous discussions on 
the problem under the guidance of Dr. C. J. Garrison. 
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II. THEORETICAL CONSIDHRATIONS 



A schematic representation of the problem under con- 
sideration is shown in Figure (1) . A rigid circular cylin- 
der of radius a submerged with center at depth d in water 
of depth h is subjected to a train of regular waves of wave 
length L and height H progressing in the positive x-direc- 
tion. The primary interest lies in the horizontal and ver- 
tical components of force induced by the waves on the 
circular cylinder. 

The complete theoretical solution to the present prob- 
lem is difficult, even assuming the flow to be unseparated; 
and, therefore, in order to proceed with an experimental 



approach, it is first necessary t 



r\ pt-v'v-kl'vr Oil'll- 



ysis to the problem. It is known a priori that the maximum 
(horizontal or vertical) force per unit length of cylinder 
is dependent upon the following variables: 



F or F = f (h ,d , a , L ,H , p , g ,y) (1) 

xmax ymax 

where h denotes the water depth, d the depth of submergence, 
a the cylinder radius, L the wave length, H the wave height, 
p the fluid density, g the acceleration of gravity, and y 
the fluid viscosity. The period of the incident wave is not 
required in the above relationship because a second rela- 
tionship already exists between the wave period and the three 
variables h, H, and L. This is true not only in the case of 
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Figure 1. Definition Sketch. 



linear waves, but in the case of large amplitude non-linear 
waves . 

A dimensional analysis of the physical quantities in- 
dicated in Equation (1) yields the following dimensionless 
variables ; 



F 

xmax 

Pga^ 



or 



F 

ymax 

Pga^ 



f 



h/a, d/a, H/2a, — p== 
L PVgha^ 



( 2 ) 



The last term on the right hand side of Equation (2) rep- 
resents the ratio of the Froude number and the Reynolds 
number and, accordingly, is an indicator of the ratio of 
viscous forces to gravity induced forces. If this number 
is large, corresponding to a large scale flow and fluids 
of small viscosity, ?t may be suspected that the gravita- 
tional forces may dominate the flow, and the viscous ef- 
fects may be considered unimportant. As a result, this 
parameter may be disregarded and Equation (2) may be writ- 
ten as : 



F F / 

xmax ymax ( 27ra , , ,/ 

or = f , h/a, d/a, H/2a 

pga^ pga^- \ L 



(3) 



A second reason for the disregarding of this parameter 
may also be cited. In unidirectional motion starting from 
rest and moving with constant acceleration, it has been 
shown by Sarpkaya and Garrison [4] that the total force 
exerted on a circular cylinder can be expressed as a func- 
tion of the relative displacement of the fluid only. Al- 
though the present problem involves oscillatory as opposed 
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f 



J 



to unidirectional motion, it is nevertheless suspected that 
the same general kind o£ dependence on relative displace- 
ment would result. In the present problem, the amplitude 
of the relative fluid motion is proportional to the wave 
height to cylinder diameter ratio and, therefore, H/2a 
represents an equivalent parameter. In the experimental 
program, an attempt is made to correlate the two force co- 
efficients on the left hand side of Equation (3) with the 
dimensionless ratios on the right hand side. 

Although the complete theoretical solution to the pro- 
posed problem is very complex, it is possible to obtain an 
approximate solution if certain simplifying assumptions are 
applied. The velocity potential associated with the regulai* 
progressive wave alone is given by: 



g^ cosh k(y+h) 

$ = —A T— -ru cos (kx-at) 

a cosh kh 



(4) 



and from this expression, the velocity and acceleration of 
the fluid particles within a avave may be written, respec- 
tively, as: 





cosh 


k 


(y+h) 


0 


cosh 


kh 


glcH 


sinh 


k 


(y+h) 


0 


CO: 


sh 


kh 


1 H 


cosh 


k 


(y+h) 


a 


CO 


sh 


kh 


vH 

gk^ 


sinh 


k 


(y+h) 


a 


CO 


sh 


kh 



sin 


(kx-ot) 


(5) 


cos 


(kx-at) 


(6) 


cos 


(kx-at) 


(7) 


sin 


(kx-at) , 


(8) 
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In wave force problems involving objects whose dimen- 
sions are small in comparison to the wave length of the in- 
cident wave, it is common to assume that the object does 
not affect the incident wave. It is then assumed that the 
object is subjected to a uniform oscillatory flow having 
acceleration equal to that induced by the \^ave at the cen- 
ter of the object if the object were not there. Following 
this line of reasoning, the horizontal force per unit length 
of cylinder may be written as the product of the virtual 
mass and fluid acceleration as 



= pTra^ 



C9) 



if it can be 
Substituting 
y=d gives the 



assumed that flow separation does not occur, 
for u froin Equation (7) and evaluating at 
dimensionless horizontal force coefficient 



as 



X 



Pga 



2 



7t(1 + C ) 




cosh ^ (h-d) 
cosh — ^ h 




( 10 ) 



Since the dimensionless force coefficient is linear in the 
wave height, this expression may be written more conveniently 
as 



X 



pga2(H/2a) 






Zira 



cosh 



2ira 



(h/a - d/a) 



cosh 



27Ta 



h/ a 



( 11 ) 



The parameter occurring in Equation (11) denotes the 
added mass coefficient and depends on the shape of the object 
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in question as well as the proximity of the bottom 
surface. Values for the configuration in question 
corresponds to a circular cylinder between a rigid 
and free surface are not available to this writer' 
ledge. The attempt in this thesis, therefore, is 
mine by trial and error the value of which give 
best fit to the experimental data. 



and free 
which 
boundary 
s know- 
to deter- 
s the 
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III. EXPERIMENTAL INVESTIGATION 



A. APPARATUS 

The wave channel as depicted in Figure (2) was de- 
signed with a rectangular cross section two feet deep with 
a width of 15 inches. The channel was fabricated in sec- 
tions. Each section was composed of two vertical walls 
and one horizontal floor member. The vertical walls were 
stiffened by use of three 4 foot long 2 x 4's attached to 
3/4 inch x 8 foot x 3 foot sheets of plywood at the ends 
and middle of its longest dimension. The channel floor 
member was composed of two 15 inch x 8 foot long plywood 
sheets separated by a 1 inch x 4 inch strip around the 
edge. The bottom and side members were mated together 
with long 3/8 inch threaded rods and formed the basic 
channel cross section. Six of these 8 foot sections were 
then bolted together to form the basic channel. Square- 
ness of the channel cross section was maintained by use 
of spacers placed between the stiffening members near the 
floor. Every butt joint was given a heavy application of 
silicone sealant to prevent leakage. 

At a location 30 feet from the wave generator end of 
the channel, a 3 foot wide plexiglass observation window 
extending the full height of the channel was installed in 
both walls of the channel as shown in Figure (3). The 
test cylinder was then mounted in the center of this win- 
dow . 
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Figure 2. Wave Channel. 




Figure 3. Wave Channel Observation Window. 
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Dissipation of the generated waves was provided by the 
beach structure located in the last section of the channel. 
The beach as shown in Figure (4) consisted of three layers 
of perforated stainless steel sheeting. Layer separation 
and structural framing v;ere provided by 2 x 4 runners with 
their largest cross section dimension perpendicular to the 
sheeting. The slope of the channel beach was maintained 
at 6:1. 



The wave generating component of the channel as shown 
in the photograph of Figure (5) had two basic elements, the 
piston-type plate and the driving motor. The plate was an 
aluminum sheet 1/8 inch thick with adjustable edges on the 
vertical sides giving it the capability of adjusting to 
conform to the sides of the channel. The plate was mounted 
on a framework attached to 4 slider bearings, and the 
slider bearings rode on two 1.0 inch brass rods attached 
to the channel walls. A driving rod ivas attached to the 
plate structural members and connected to a flywheel at- 
tached to the output shaft of a Vari-Drive motor. The 
drive rod was attached to the flywheel by means of a bolt 
and slot in such a manner that the eccentricity could be 
varied between zero and 6 inches. In this way, waves of 
various amplitudes could be generated. The variable speed 
transmission was driven by a 2 hp inductance motor and pro- 
vided a speed range extending from 22 to 180 rpm giving 
wave lengths ranging from 18 feet to 1.5 feet at a water 
depth of 18 inches. 
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Figure 4. Views o£ Channel Beach. 
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Figure 5 . Wave Generator 
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A plexiglass circular cylinder o£ 6 inch diameter with 
an end clearance of approximately 1/16 inch was used in the 
test reported herein. As depicted in Figure (6) , the cylin- 
der was hung by two fine wires at the appropriate distance 
off the wave channel floor. Weights were added inside the 
cylinder to pretension the ^^^ires and to prevent slack in 
the wires during testing. 

The upper ends of the wires were attached to small 
aluminum cantilever beams equipped with strain gages on 
both the compression and tension side of the beams. These 
vertical force strain gages were then connected into a 
Wheatstone bridge and read-out on one channel of a two- 
channel Sanborn Model 296 amplifier/recorder. The cylinder 
was restrained in the horizontal direction by similar fine 
wires connected into a harness and around a 4 inch diameter 
ball bearing pulley mounted on the v/ave channel floor as 
shown in Figure (7) . These wires were then attached to a 
pair of cantilever beams equipped with strain gages, and 
the strain gages were connected into a Wheatstone bridge 
and read-out on the second channel of the two-channel am- 
plifier/recorder. 

The four beams, two for the vertical force and two for 
the horizontal force, were attached to a plywood mounting 
board positioned over the channel as shown in Figure (8) . 

The beams were all exactly alike, and the section of the 
beams where the strain gages were mounted was reduced to a 
cross section of 1/4 x 1/4 inches from a nominal beam cross 
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Figure 6. Cylinder Suspension System. 




Figure 7. Horizontal Wire Harness and Pulley Assembly. 
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Figure 8. Beam Mounting Board. 
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section o£ 1/2 x 3/8 inches. The machining provided a 
smooth transition into the area of stress concentration. 
SR-4 FAE-25-12S13 strain gages were applied to the top 
and bottom surfaces of the beams. Each beam was 9 inches 
long and had mounting holes for two moment arm positions, 
one at 6 inches and one at 3 inches from the strain gage 
location. This variable load attachment point gave each 
beam a sensitive force measurement capability under a 
large range of values. The beams were designed to with- 
stand a six-pound load at the largest moment arm. 

The wave height gage, as shown in Figure (9), con- 
sisted of a small diameter aluminum rod acting as the 
basic frame to which a plexiglass mounting table was at- 
tached. The elements of the Wheatstone bridge were 
mounted on this table, and two 30-gage copper wires spaced 
1/2 inch apart were extended from the table and attached 
to an insulator mounted at the end of the aluminum rod. 

An 18-inch separation between the levels of the table and 
insulator provides a large measurement range capability. 

The other end of the rod was mated to a traverse block 
stand which allowed the gage to travel ± 6 inches, and 
provided a convenient, yet accurate, means for calibration 
by adjusting the submergence of the wires. The complete 
unit was positioned over the channel and centered at ap- 
proximately 5 feet in front of the cylinder. 

The wave height probe was of the parallel wire re- 
sistance type. This type of gage operates on the principle 
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Figure 9. Wave Height Gage. 
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that the conductance between the two parallel submerged 
wires varies approximately proportionately to the length 
of wire submerged. The wave gage was connected in parallel 
to one leg of a Wheatstone bridge as depicted in Figure 
(10). To compensate for the inability of the carrier am- 
plifier to offset the resistance of the immersed parallel 
wires, an external variable lOK ohm potentiometer was con- 
nected across an opposite leg of the bridge. This resistor 
was used for coarse balancing of the bridge. The Wheatstone 
bridge was connected vip signal and excitation leads to a 
portable carrier preamplifier and read out on a one-channel 
H/P Sanborn 301 amplifier/recorder. 

During preliminary operation while attempting to clear 
the problems from the system, it was found that installation 
of an impermeable wall between cylinder and channel floor 
was a necessity. The experiment was designed to simulate 
the placement of the circular cylinder on the bottom, but 
it was necessary to maintain a small but finite gap between 
cylinder and the bed of the channel to insure that the 
pair of wires connected to the cantilever beams provided 
the sole vertical support. However, it was found that the 
vertical forces as well as, to a lesser extent, the hori- 
zontal forces were strongly dependent on the gap width. 

Flow visualization studies showed that a high velocity 
jet occurred in the small gap under the cylinder. Thus, 
in order to eliminate the flow under the cylinder and yet 
not influence the force measurement system, a flexible 
plastic slieet was installed across the gap between the 
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Figure 10. Wave Height Gage Bridge Scliematic. 
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cylinder and floor. The sheet was joined to the cylinder 
by pressing the plastic into a machined slot in the cylin- 
der wall, and then anchoring it with a length of "0"-ring. 
The other end of the plastic was laid under a metal sheet 
12 inches x 14 1/2 inches and 1/16 inch thick and placed 
at the channel bottom as shown in Figure (11) . The gap 
between the cylinder and floor was adjusted to approxi- 
mately 1/8 inch. This gap was found to be large enough 
to insure that the cylinder never touched the bottom when 
subjected to the largest wave. 

In order to demonstrate the rather large influence of 
this small gap under the cylinder. Figure (12) was prepared 
sho^^?ing both the maximum horizontal and vertical force as 
a function of wave amplitude. It is apparent from this 
figure that removal of the plastic barrier had the effect 
of reducing the maximum horizontal force a slight amount. 
With the barrier in place, the upward force was generally 
quite large while the downward force was rather small. 
However, upon removal of the barrier, a rather high veloc- 
ity flow through the gap was observed by injection of dye 
near the wave channel floor, and this local region of high 
velocity apparently resulted in low pressures on the under- 
side of the cylinder causing the downward force to be in- 
creased considerably. The rather extreme effect of the 
gap is clearly demonstrated by comparison of the various 
curves indicated on Figure (12) . 
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Figure 11. Barrier Assembly and Wave Profile. 
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Relative Wave Height, H/2a 

Figure 12. Comparison o£ Wave Force Coefficients Shoving 
the Effect of the Barrier. 
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Inasmuch as the intent of the present study was to 
determine wave forces on a cylinder lying on the bottom 
wherein the gap would not be present, the plastic barrier 
was maintained in place throughout the test program for 
all cases where the cylinder was placed near the bottom. 

B. TEST PROCEDURE 

The first step in the test procedure was to fill the 
wave channel with water to obtain the desired water depth 
to cylinder diameter ratio, h/a. Experiments were con- 
ducted at the three depths, 12, 15, and 18 inches. After 
the desired water level was attained, a careful inspection 
was carried out to insure that all equipments were properly 
prepared for operation. 

The vertical wire position scre^^?s were adjusted to 
obtain the desired elevation of the cylinder off the chan- 
nel bottom. The horizontal wires were pretensioned to 
prevent slacking and the attendant errors in force measure- 
ments . 

A typical run may be described as follows. With water 
adjusted to the proper depth, a given wave length, or equiv- 
alently, wave period was established by setting the speed 
of the Vari-Drive unit. Subsequent to a sufficient warm-up 
period, the amplifier/recorders shown in Figure (13) were 
balanced and readied for calibration. Calibration of the 
three recording channels, horizontal force, vertical force, 
and wave height was carried out at the attenuation expected 
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Figure 13. Experimental Equipment. 




to be in use during the run. After calibration, the wave 
height was set by adjusting the eccentricity of the wave 
generator to the appropriate amplitude. The wave generator 
was then turned on, and after a steady state wave was es- 
tablished in the channel, the recorders were turned on and 
data was taken. The wave generator was then switched off, 
and the next wave height was set on the flywheel eccentric 
while maintaining the same speed setting. Wave heights at 
the fixed wave length were varied within the capability of 
the system to produce representative readings usually ranging 
from approximately 1/4 inch to 5 inches. At the conclusion 
of each series of runs at a given wave length, a check cali- 
.bration on the three channels was quickly carried out. Then, 
a new wave length was selected and set on tlie Vari-Brive, 
and the cycle of balancing, calibration, wave height selec- 
tion, and wave production was repeated. 

Calibration was performed on the load cells as follows. 
The recording chart was set at I'^'^/sec. speed, and loads in 
one-pound steps were placed on the saddles located beneath 
and hung from the horizontal sensing wires as shown in Fig- 
ure (14). Calibration loads were placed initially on the 
after beam, then on the forward beam. The vertical load 
cell was calibrated by use of a platform of known w’eight 
centered on the cylinder as shown in Figure (15) . Weights 
were then placed on this platform. 

Wave height gage calibration was performed by adjusting 
the vertical position of the gage through the use of a slider 
bearing and screw. The procedure was carried out as follows. 
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Figure 14. Horizontal Force Calibration Operation. 
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Figure 15. Vertical Force Calibration Operation. 
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The gage was raised 1/2 inch, simulating a decrease in water 
level, then lowered 1/2 inch below the null, simulating an 
increase in water level. Then, the gage was raised +1 inch 
and lowered -1 inch, etc. This procedure closely simulated 
the actual passage of a wave profile with crests and troughs. 
The generated wave crest passing through the test location 
is shown in Figure (11) . 
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IV. EXPERIMENTAL RESULTS 



As indicated in Equation (3) in the THEORETICAL CON- 
SIDERATIONS, the wave force coefficients can be represented 
as functions of relative water depth, relative depth of sub- 
mergence, relative wave length, and relative wave height pa- 
rameters. Accordingly, a test program was devised to measure 
the wave force coefficients for practical ranges of the in- 
dependent parameters. Four series of tests were conducted 
at different water depths and depths of submergence. Three 
of these corresponding to depth to cylinder radius ratios 
of 4.0, 5.0, and 6.0, were carried out with the cylinder 
near the bottom, i. .e., at d/a= 3.0, 4.0, and 5.0, respec- 
tivcly. A four til sorios of tosts conJuctod with tho 

cylinder raised one radii off the bottom corresponding to 
d/ a= 4.0 at a water depth ratio of 6.0. 



A. TYPICAL WAVE FORCE TRACES 



Inasmuch as it is difficult to completely characterize 
the force variations with a single number, typical strip 
chart read-outs for the various series of tests are pre- 
sented in Figures (16-28). Figures (16-19) show the hori- 
zontal and vertical force variations with the cylinder on 
the bottom corresponding to a relative depth of submergence 
of d/ a = 5.0 and a relative water depth ratio of h/a= 6.0. 
The figures are arranged according to increasing values of 
the wave length parameter ranging from 27ra/L= 0.09 to 0.97. 
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On each figure, data corresponding to four values of rela- 
tive wave height are presented in ascending order. 

Figure (16) representing the largest wave length and 
having a depth to wave length ratio of h/L= 0.08 corre- 
sponds to waves which are somewhat outside the sinusoidal 
range. These long waves are characterized by long, smooth 
troughs and sharp crests. Accordingly, the horizontal 
force traces are characterized by positive and negative 
peaks which result from the maximum acceleration occurring 
on either side of the crest. These two peaks are followed 
by a relatively flat portion of the force trace corres- 
ponding to the long, smooth trough portion of the wave. As 
can be seen from the figure, this kind of variation appears 
to predominate even as the wave height is increased foi' 
this long wave case. 

The vertical force, however, shows a much more incon- 
sistent variation. If the system natural frequency is dis- 
regarded, definite upward peaks are observed at points on 
the trace corresponding to the wave crests. These peaks 
are apparently related to the reduced pressure over the 
upper portion of the cylinder ' which is associated with the 
velocity squared term in the Bernoulli equation since the 
maximum velocity occurs at the crests. In sinusoidal waves, 
an equal and opposite velocity occurs at the troughs which 
would naturally give rise to a second upward peak in the 
force trace. However, in shallow water waves, the veloc- 
ity associated with a trough is much reduced and, there- 
fore peaks corresponding to the crests only are observed. 
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Figure (17) corresponds to a shorter wave length having 
a relative wave length value o£ 2ira/L= 0.16. It is apparent 
from the figure that both the horizontal and vertical forces 
vary sinusoidally for the smallest amplitude wave correspond 
ing to H/2a= 0.5. As the wave height is increased, it is ap 
parent that the horizontal force continues to vary sinusoid- 
ally for the full range of amplitudes. The vertical force, 
however, shows a sinusoidal variation for the smallest am- 
plitude waves and becomes increasingly less sinusoidal as 
the \\^ave height is increased. Moreover, two separate peaks 
appear to develop on the vertical force trace as wave height 
increases, one associated with the maximum velocity occur- 
ring at the crest and the other associated with the trough. 
The maximum velocity at botli the crest and the trough re- 
sults in reduced pressures on the upper portion of the cylin 
der and, consequently, corresponds to upward vertical forces 
Figures (18) and (19) correspond to still smaller \vave 
lengths of 2ira/L= 0.24 and 0.77, respectively. It may be 
noted that again the horizontal force shov\fS a sinusoidal 
variation throughout the range of wave amplitudes while the 
vertical force sho^vs the sinusoidal variation at small am- 
plitudes becoming non-sinusoidal as the amplitude is in- 
creased. Figure (19), however, which corresponds to the 
shortest wave lengths shows sinusoidal variations through- 
out the range of wave amplitudes for the horizontal as well 
as vertical force. 
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Horizontal and vertical force coefficients correspond- 
ing to relative water depth of h/a= 5.0 with the cylinder 
on the bottom, i.e., d/a= 4.0 are shown in Figures (20-22). 
These three figures represent three different values of 
relative wave length extending from 2Tra/L= 0.10 to 0.67. 

The general trends tend to be similar to those corresponding 
to the larger depth of h/a as shown in Figures (16-19) and 
previously discussed. Data corresponding to the smallest 
value of 2iTa/L as presented in Figure (20) shows trends 
typical of nonlinear shallow water \;aves . The horizontal 
force shows the null period associated with the long, smooth 
troughs, while the vertical force is characterized by the 
upward force peak corresponding to the wave crests. As the 
v.’ave length is decreased, the horizontal force again dis- 
plays sinusoidal variation as indicated in Figures (21) and 
(22). For the shorter wave lengths, the vertical force shows 
a sinusoidal variation at small amplitudes with superimposed 
upward peaks occurring near the maximum velocity location at 
the crests and troughs. Figure (22) corresponding to the 
shortest wave length shows a sinusoidal variation for both 
force components throughout the range of wave heights. 

Typical force variation records for the series of tests 
conducted at h/a= 4.0 are presented in Figures (23-25). 

Three representative values of relative wave length, 2iTa/L= 
0.11, 0.21, and .45 are presented. Figure (23) shows varia- 
tions similar to those of equal values of 27ra/L at the two 
larger depths already discussed. Figures (24) and (25) arc 
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also similar to previous runs, the primary difference being 
the increased magnitude occurring at the smaller depths. 

The final series of runs was carried out at a relative 
water depth ratio of h/a= 6.0 with the cylinder suspended 
one radii off the bottom, i.e., d/a= 4.0. Typical force 
traces at three different values of relative wave length, 
27ra/L= 0.10, 0.16, and 0.6, are presented in Figures (26), 
(27), and (28), respectively. The horizontal force trace 
as shown on the figures reflects the same general trends as 
the previously described series. The vertical component 
of force a] so shows trends similar to those previously dis- 
cussed with the cylinder on the bottom, but the peaks as- 
•sociated with maximum velocity in the wave occurring at the 
crests and trouglis are seen to be in the opposite direction. 
That is, the velocity of the fluid through the one radii gap 
between the cylinder and the floor is apparently consider- 
ably larger than the velocity through the space above the 
cylinder and beneath the free surface and, therefore, a 
sucking-down action as opposed to a lifting action results. 
As in previous runs, it may be noted that Figure (28), cor- 
responding to the largest value of 2va/L, shows a sinusoidal 
variation for both force components throughout the range of 
wave heights. 

B. WAVE FORCE COEFFICIENTS 

In order to represent the results of the experiments in 
a simpler, albeit less descriptive form, the maximum values 
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of the dimensionless force coefficients obtained from the 
force traces were plotted as a function of relative wave 
height, H/2a. Results corresponding to the relative water 
depth of h/ a= 6.0 with the cylinder on the bottom, i.e., 
d/ a= 5.0, are presented in Figures (29-32). Results cor- 
responding to the other two relative water depths of h/ a= 

5.0 and 4.0 also with the cylinder on the bottom are pre- 
sented in Figures (33-36) and (37-40), respectively. These 
results correspond to values of 2iTa/L ranging from approxi- 
mately 2ira/L= 0.08 to 0.97. 

As was apparent from the force traces, the horizontal 
force showed sinusoidal variations with equal magnitude in 
both directions and, therefore, only one curve was required 
to represent the amplitude of this variation with relative 
wave height. The dimensionless force coefficients are de- 
fined as f=F/pga^, and three curves are noted on each fig- 
ure, one for the horizontal and one each for the upward and 
downward vertical force. As is evident from the figures, 
the horizontal force coefficient shows a linear variation 
with the wave height in all cases. 

However, as the vertical force showed, in general, larger 
upward than downward magnitudes, the two curves labeled f^ 
and fy are presented showing the maximum upward and downward 
values of the vertical force coefficient, respectively. 
Throughout the range of the data, the vertical force coeffi- 
cient for the cylinder located on the bottom showed a quite 
consistent variation, the upward force being always larger 
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Figure 29. Wave Force Coefficients for h/a= 6.0, d/a= 5.0. 
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Figure 30. Wave Force Coefficients for h/a= 6.0, d/a- 5.0. 
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RELATIVE WAVE HEIGHT, H/2a 

Figure 31. Wave Force Coefficients for h/a= 6.0, cl/a= 5.0. 
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RELATIVE WAVE HEIGHT, H/2a 

Figure 32. Wave Force Coefficients for h/a= 6.0, d/a= 5.0. 
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than the downward force and having a generally concave up- 
ward trend. The downward force was always quite small show- 
ing a nearly linearly increasing variation followed by a 
breaking off at higher wave amplitudes. 

The variations of the force coefficients as observed 
in Figures (29-40) can be justified if one considers the 
components which make up the dynamic pressure acting on the 
surface of the cylinder. Two components of pressure are 
generally recognized, the one component being associated 
with the unsteady term in the Bernoulli equation, and the 
second being associated with the velocity squared term oc- 
curring therein. The unsteady term involves tlie velocity 
■potential to the first power and, therefore, is linear in 
H/2a. The second term involves llie spatial derivative of 
the velocity potential squared and is, therefore, propor- 
tional to (H/2a)^. Thus, linear variation of the horizon- 
tal force with H/2a apparently results from contributions 
from the unsteady component of pressure only. The velocity 
squared term which would ^contribute a nonlinear component 
of force can contribute nothing to the horizontal force if 
the flow pattern is symmetrical about the y-axis. This, 
of course, is expected provided separation does not occur 
and accounts for the linear variation in the horizontal 
force with the parameter H/2a. 

The vertical force variation with H/2a, however, shows 
markedly different trends. Both the upward and downward 
components seem to vary linearly with relative wave heicht 
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Figure 33. Ivav'^e Force Coefficients for h/a= 5.0, d/a= 4.0. 
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RELATIVE WAVE HEIGHT, H/2a 

Figure 34. Wave Force Coefficients for h/a= 5.0, d/ a= 4.0 
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Figure 35. Wave Force Coefficients for h/a= 5.0, d/ = 4.0. 
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RELATIVE V/AVE HEIGHT, H/2a 
Figure 36. Wave Force Coefficients. 
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RELATIVE V/AVE HEIGHT, H/2a 

Figui'e 37. Wave Force Coefficients for h/a= 4.0, d/a= 3.0. 
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RELATIVE WAVE HEIGHT, 11/ 2 a 

Figure 38. Wave Force Coefficients for h/a= 4.0, d/; = 3.C. 
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RELATIVE WAVE HEIGHT, H/2a 

Figure 39. Wave Force Coefficients for h/a= 4.0, cl/a= S.fi. 
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Figure 40. Wave Force Coefficients for h/a= 4.0, d/a- 3.0. 
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for small amplitude waves, and the upward force tends to 
bend upwards while the downward force tends to decrease at 
larger values of H/2a. These variations, however, are con- 
sistent with the previous discussion regarding the pressure 
on the surface of the cylinder. In the case of the vertical 
force, the flow pattern is not symmetrical about a horizon- 
tal axis passing through the center of the cylinder, and, 
therefore, the contribution from the velocity squared term 
in the Bernoulli equation does not cancel. As a crest or 
trough passes the cylinder, the wave induced velocity is a 
maximum giving rise to a component of pressure (or force) 
proportional to (H/2a)^. The vertical force coefficient 
•variation with H/2a, therefore, is composed of two terms, 
one proportional to 11/ 2a and one proportional to (H/2a)^. 
Thus, its somewhat parabolic shape is not unexpected. 

The maximum values of the horizontal and vertical force 
coefficients for the final series of runs conducted at a 
water depth ratio of h/a= 6.0 with the cylinder suspended 
off the bottom, i.e., d/a= 4.0, are presented in Figures 
(41) through (44)., As in the previously described tests, 
the horizontal force coefficients showed a linear variation 
with wave amplitude whereas the vertical force coefficients 
showed an initial linear portion followed by a nonlinear por- 
tion. However, in this series of runs, the high velocity 
in the gap beneath the cylinder caused a net nonlinear con- 
tribution to the force in the downward direction making the 
magnitude of f^, larger than f , a trend opposite to tiio>c 
associated with the cylinder near the bottom. 
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Figure 41. Vlave Force Coefficients for h/a= 6.0, d/a= ■!.'). 
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Figure 42. Wave Force Coefficients for h/ a= 6.0, d/a= 4.0. 



79 



if 

% 



1 



FORCE COEFFICIENTS 





Figure 43. IVave Force Coefficients for h/a= 6.0, d/a 4 . J . 
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Figure 44. Wave Force Coefficients for h/a= 6.0, d/a= 4.0. 
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HORIZONTAL FORCE COEFFICIENT 



The maximum values o£ the horizontal force coefficient 
f , as shown in Figures (29-40) , are seen to vary linearly 
with the relative wave height, H/2a. In view of this lin- 
ear behavior, it is possible to represent these forces in 
a much more condensed form. That is, since the force co- 
efficients are linear with H/2a, the parameter. 



F^/[pga2(H/2a)] , 



which actually represents the slopes of these curves can 
be presented in a single figure as a function of the rela- 
tive wave length, 2ira/L. Such a plot is presented in Fig- 
ure (45) wherein this force coefficient is represented as 
a function of 2iia/L for tlie four senes of tests earned 
out at different Avater depths and depths of submergence. 

The upper three curves shoA^n in Figure (45) represent 
three different ANrater depths Avherein the cylinder Avas lo- 
cated on the bottom. The effect of A\rater depth is clearly 
indicated, as the curves corresponding to decreasing depths 
lie above the one at a maximum A'/ater depth of h/a= 6.0. 

The general trends of these three curves are all similar, 
maximum values of force coefficients occurring near 2Tra/L= 
0.3 A^^ith decreasing force at shorter Avrave lengths. At 
longer A^:ave lengths than those corresponding to the maximum 
value for the force coefficient, the curves tend to become 
flat Adiile the linear theory A\rould indicate a decrease 
similar to that associated Avith the curve corresponding tj 



82 



FORCE COEFFICIENT, F / [ pga^ (H/2a) ] 




Relative Wave Length, Zua/L 
Figure 45. Horizontal Force Coefficients. 
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h/ a= 6.0 and d/a= 4.0. It is possible that this behavior 
is caused by the fact that the wave is becoming a nonlinear 
shallovi? water wave. 



The horizontal force coefficients corresponding to the 
fourth series of tests carried out at a water depth ratio 
of h/a= 6.0 and with the cylinder suspended at one radii 
off the bottom is also shov>:n in Figure (45) as the lowest 
curve. The general trends are seen to bo somewhat similar 
to the cases where the cylinder was located on the bottom 
except that the curve tends to be somewhat flatter. 

In an attempt to compare the experimental results with 
the approximate theory, Equation (11) was programmed on a 
digital computer, and values for the horizontal force co- 
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various values for the added mass coefficient, C . These 
results were then plotted and compared with the experimental 
results in order to determine which value of gave the 
best fit. These curves representing the "best fit" were 
then plotted on Figure (45) and were compared v/ith the ex- 
perimental plots. As might be expected, it was possible 
to obtain good agreement between the approximate theory 
and the experimental data for the largest depth ratios. In 
these cases, the proximity of the free surface, which was 
disregarded in the approximate analysis, was probably not 
particularly important. However, for the shallowest depth 
where the free surface was more influential, it appeared 
to be impossible to find a value for vvrhich gave reason- 
able agreement with the plotted experimental results. 
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One point of interest in Figure (45) was to compare 
the two curves associated with a maximum depth of h/a= 6.0, 
the one case corresponding to the cylinder located on the 
bottom, d/ a= 5.0, and the other witli the cylinder elevated 
one radii off the bottom, d/a= 4.0. As can be seen from 
the figure, the latter curve is considerably below the for- 
mer and corresponds to an added mass coefficient of approxi- 
mately C^= 1.2. The data corresponding to a similar 
situation with the cylinder located on the bottom appears 
to agree relatively well witli the curve based on = 2.8. 
This marked increase in corresponding to the case where 
the cylinder was located on the bottom apparently results 
from the proximity of the bottom. 

A second inteiesting feature in Figure (45) was the 
apparent discontinuity in the two curves associated with 
the two shallowest steps near 27ra/L= 0.3. The cause of 
this discontinuity is not certain, and it might be dis- 
missed as scatter in the experimental data. Tliere is, 
however, a second explanation which is equally probable. 

It is possible to have a small amount of reflection at 
this value of 2Tra/L, and it is possible that this reflec- 
ted waves affected the wave height measurement. Depending 
on the phase of the reflected wave, it may eitlier add or 
subtract and cause an apparent increase or decrease in the 
value of the wave force coefficient. The absence of tliis 
discontinuity in the curves corresponding to the greater 
depths is in agreement with this liypothesis since reflec- 
tion decreases with relative deptli. 
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V. CONCLUSIONS 



An experimental study of gravity wave interaction with 
a submerged horizontal cylinder has been conducted. From 
the results of this study, the following conclusions are 
drawn . 



1. For the range of wave. heights tested in the experi 
ment , the horizontal forces resulting from gravity wave 
interaction with the cylinder increased linearly with wave 
height . 



2. Tlie vertical force was alv.’ays much smaller than 
the horizontal force when the cylinder was located on the 
bottom. Tlie downward part of this force fluctuation was 

crease witli \\'ave amplitude. At very small amplitudes 
where the velocity squared contribution to the pressure 
was small, a sinusoidal variation of the vertical force 
resulted . 



3. For the case corresponding to the cylinder sus- 
pended one radii off the bottom, the horizontal force was 
reduced and the downward force was greater than the upward 
force . 



4. For cases where the depth of submergence was large 
the approximate analysis appeared to give valid results 
provided the proper value of the added mass coefficient was 
hnown . 
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5. When the value of Zva/L becomes small, correspond 



ing to long waves 
similar to those 
In particular, a 
occurs during the 



, the wave force trace shows features dis 
associated with the more sinusoidal waves 
null period in the horizontal wave trace 
passage of the trough. 
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TABLE I - continued 
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4 h 
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*H •••••9990990999999999999 

^ro m CO CM ^ —I O o CO CM rH O o CM CM t M CM CM CM ^ OM rM r-* 3 
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TABLE II - Experimental Wave and Force Data, h/a= 5.0, d/a= 
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00 CM or^ r-f m nO r-H 00 00 o CM o r-H nO 00 m 4- 00 *-« o CM o nO 4- nO cn CM o CM m r-j m 
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[_) 0 NOsONO'ONOxONO>OQOcor^r^oOr^r^r^4‘4’4'4'4*4‘4’CMcOCM(McMCMCMCM30v.^OO 
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TABLE II - continued 
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TABLE III. Experimental Wave and Force Data, h/a= 4.0, d/a= 
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yf —KM^vOirv^cMr-HOcor^ ^om'4■m^-lC^cD^-NOlr^M■cM.~lC^oo^o^A^4■^ocMoa^oo^o^A 
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oooooooooooooooooooooooooooooooooooo 



, c>cM'OLnh-cocMsOCDtnocMcocMuoLncoo-HvOCMOcM^cM>tocOrHr--ocr'inrococ> 
^ c^co^^co^-lncocMv^{^^cMOC^^-cOr-l^O'ncMoo^-m^O(M^-^c^^-vnl^cMO^•(^^oo 
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V y 1—4 I-H r-H C? C5 (O O* CO CD CiJ 'O iO CD CD O 

^ '^*>^o^oini-Hcoc^o^'OsOinoocMo>r-HM-o o^-NtrHoc^'4■a'sO'OLn^4•vO^.D'tlA 
coinr^r-Hinocoo^M“(7^vOCMcocor^CMnror^r-HOOroO'4*!X)OLO -4U^c^'?■ln(M'^ao 

IT ■ . «®00*#0oe#«0e#0oe00«#0#o0og**e«^ce®0 

inxtcOrOCMCMr-HO CtOCOcOCMCMr-HO"^’>^'COcOCM.-Hr-HON5'fOrOCMCM-HOONf M-cOCM 
^ J LnrHcocxDC^ooLn^o-^cMCMr^•-4sOOr-HvY^lncoc^^-lCMno^^covOv^ OLn^o^cxDCM^'0 

^-stin.-HLO^icxDvO C^or^'4-coc^>i’r^oaDino<^cMoouo-Tcooo<f‘cocCin(M— icocvjin 

♦•♦•••e«0eg««g«ee«0g0«»e«o*«««*4»6oo« 

cOCOCMCM'-Hr-HOOCOcOCMCMCM-H.-HOcorMcMCM'-HrHO O CM CM OOO^tcOcOCM 

r, ^ in^oh-coaDooi^cOr-HO‘-H'^'-<'~<'-«Of^mvO '0'OsONCr^o^oa*ocooocoooN^-M'>^-'^ 

^ vO O O sO 'O nO sO O' O O" (7^ CM CM CM CM CM CM CM CM r-H .-H.-H 

CMCMCM CM CMCM CM (M»-Hr-H.-H,-H'-HrHr-Hr-4.-Hr-H.-Hr-Hi-Hr-lf-Hr-Hr-H^i-Hr-Hr-Hr-H»-H.-H.-Hr-lr-4.-H 



co^tm sO h-co(j^ooo^.Mco vMT. or- o O*^ CM TOM’ msOr^cCG^ o—i CM CO or- CO 

^ oo<x)coc30cxDcocooa'a>cro^oc>c>cr CT'OO OOOOOO 
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TABLE III - continued 



O ro <) CT' (\j in in ro c\i i-H fH rn CO 00 f\j in rH (\j r-H 00 fT) in c\j (\j fo f\i 

^ OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 



CO 00 tn (Nj in CN 'O iTi m r -4 00 00 o o m IN 00 00 00 f'- ro ovj >o 

^ ^ i-^OOrnfOCNCNr-HOrnrnfnfnC\J(\Jr-t»-H(\JC\J(\J(\JC\Jf\Jr-Hr-4»-Hr--lr-lr-l»-HO 
^ OOOOOOOOOOOOOOOOOOOCJOOOOOOOOOOO 



o r- > 1 * 00 o m m m O' 00 o oo vO vO m >t m fo nj O' r- rg a> 00 

>t^n-^^-^Olr^lnfnr-Hoo^-^-^o^>tfOr--lN^''4•>^*^^‘N^fnc^J^nfn(\l^^Jr^Jr-^o 

. «0cc«««>«9«>0 0*9«000009«ft0«oe0«0«« 

OOOOOO-^OO DO OOO OOOOOOOOOOOOOOOOO 
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TABLE IV. Experimental Wave and Force Data, h/ a= 6.0, d/a= 
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TABLE IV - continued 
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TABLE V. Wave Forces With and Without Barrier, h/a= 6.0, d/a= 
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WAVE FORCES ON A SUBMERGED HORIZONTAL CYLINDER" 
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